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Abstract 

CN , Experimental observations suggest that there are differences between the behavior of particles 

' jy^" opticaUy trapped in air and trapped in a liquid phase. We present a modified version of Mie 

o ■ 

Debye Spherical Aberration theory to numerically simulate such optical system in attempt to 

^: 

O , explain and predict these effects. The model incorporates Mie scattering and focussing of the 

O ' trapping beam through media of stratified refractive index. Our results show a geometrical optics 

J>-j' approach cannot correctly describe our system and that spherical aberration must be included. 

Qh. We successfully qualitatively explain the observed phenomena and those of other authors, before 
discussing the limits of our experimental techniques and methods to improve it. We draw the 

> 

' important conclusion that when optically trapping aerosols the system does not behave as a true 

o ■ 

. 'optical tweezers', varying between levitation and single beam gradient force trapping depending 



o 
o 



X 



on particle and beam parameters. 
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I. INTRODUCTION 



Optical manipulation has matured considerably in recent years to become a powerful 



quantitative tool 



m 



The forces imparted are ideally situated to enable them to act as 
force transducers in molecular biology |3|. Also their ability to trap and isolate single or mul- 
tiple objects makes them powerful in many different disciplines including biochemistry jj], 
fiuidics 5|, condensed matter and fundamental physics p, l7|. One new application area, 
which this paper discusses, is the use of optical traps for studying aerosols and their basic 
physical and chemical properties js]. 

In all these disciplines it is useful to have an understanding of the processes occurring 
including the optically generated forces, which clearly underpin many of the measurments 
being undertaken. Modelling of optical forces is already used extensively in the field of optical 
manipulation to understand, for example, force mapping ^] and optical binding [lO| or to 



extract physical parameters not otherwise obtainable from measurements 



One of the 



best examples is the understanding of how optical forces affect the cytoskeleton. Changes in 



the elasticity of this matrix are directly related to stages of cancer within individual cells [12 1. 
As we will see, when optically trapping aerosols, we are working at the very edge of what 
can be considered to be 'optical tweezers' and so the modelling of the forces involved may 
act as a method for testing the theories at their limits. 

Aerosols are a significant constituent of the atmosphere and a major factor in determining 
its chemical balance, for example the ozone hole and acid rain 13||, impairing visibility and 
contributing to radiative balance 1^. Furthermore, understanding aerosol behaviour, how 
they enter and interact with the body is relevant for both the effects of pollution on the 
human population and the efficacy of medicinal drugs 15|. 



Aerosols can be grouped into three main categories or modes. The nucleation mode 
consists of small emitted or newly nucleated particles with a mean radius less than 0.05 fim. 
Upon growth and coagulation these particles move into the accumulation mode with radii 
between 0.05 fim and 1 um. Finally the coarse mode of aerosols consists of particles with 
radii greater than 1 fim 14 1. 

In terms of surface area to volume ratio, accumulation mode aerosol constitutes the 
largest proportion of atmospheric aerosol and therefore dominates atmospheric aerosol chem- 
istry [13i]. This mode is also able to penetrate deep into the lungs playing a major role in 
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the interaction of aerosol with the human body. Its near optical wavelength size also affects 
visibility [l^. Most optical manipulation work to date concentrates on coarse mode liquid 
phase aerosols, particularly relevant in drug delivery and atmospheric chemistry. However, 
the questions that can be answered with current techniques are limited so they will need to 
be altered to allow accumulation mode aerosol to be studied. 

Our study aims to investigate the optical forces imparted to optically trapped airborne 
microscopic spheres. This is a very difficult task to carry out experimentally [igI and the 
standard technique is to computationally model momentum transfer from focussed beam to 
particle. 

Our aim in developing models examining the optical forces involved in airborne optical 
traps is to increase the understanding of the outcomes from relevant experiments to give 
a more complete picture of the process. We can also probe the boundaries of the current 
experiment to see if they can be extended. 

Symmetry within the system simplifies the mathematics. A sphere in an axially symmet- 
ric beam is probably the simplest of formulations, with a large amount of the constituent 
work already available. Should one like to model the trapping of non-spherical objects then 
the computation becomes more complex with the T-matrix approach becoming the favoured 
method. Once the T-matrix has been calculated for a given object it need not be calculated 
again for every orientation of the object in the beam, thus making it rather advantageous. 
Nieminen et al. [l^ have used this approach to code an 'optical tweezers computational 
toolbox' freely available for use |l8 |. 

The principle problems with many of the approaches available are the over-complexity 
(ours is a relatively simple problem) , their inappropriateness for the size scale we are looking 
at, and the lack of a description of the true trapping beam profile, as we will discuss shortly. 
A lucid introduction to the inadequacies is given by Viana et al. [19]. The microdroplets we 
are studying are ^ 1 /xm in radius so the force calculation lies above the Rayleigh regime [20 1 
and to a first approximation can be described by a geometrical optics (GO) model 2l|]. It 
will quickly become apparent that by studying airborne objects a GO model cannot deal 
with the system under study. To resolve this problem we implement a model that uses an 
integral representation of focussed light crossing refractive index interfaces and an exact 
form of plane wave scattering from spheres. 

We have previously shown [22] that the motion of optically trapped droplets can be 
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described using a simple harmonic oscillator model. Evidence also suggests that the odd 
phenomena observed in the optical trapping of aerosols originates from the water droplets' 
interaction with the optical field. The principle aim of this investigation is to see if the 
isolated physics of optical forces leads to the phenomena observed in experiments and to 
see how far the boundaries of the techniques can be pushed. We will begin by stating the 
experimental phenomena observed by ourselves and others: 

1. As trapping laser power increases so does the height above the water layer or coverslip 



that the droplet is trapped 



23|. 



2. With further increases of power the droplet is lost from the trap. This does not always 
occur and is more pronounced for smaller droplets 22|. 



3. After first capture, the droplet can undergo significant growth or evaporation coupled 



with large axial oscil 



but are far slower 



24 



ations. These oscillations can occur significantly after capture 



25|. 



4. There is a linear dependence of 'captured' dro 
droplets cannot be trapped with high power [26 



et radius with trapping power. Small 
2'; 



The boundaries we would like to determine and push are; 

• the range of particle refractive indices it is possible to trap, 

• the limits of particle radius that can be trapped with the current apparatus. 

Fortunately our system is highly symmetric and almost complete rigorous wave theory 
solutions already exist in literature, although no computer code is readily available and 
modifications must be made to suit our problem. 

Here the forces exerted on spheres are decomposed into two directions, lateral and axial. 
That is the direction perpendicular to and crossing the beam propagation axis, and the 
direction lying on the beam propagation axis respectively. The simulations programmed in 
MATLAB calculate, for a given point on one of the two axes described, the efficiency with 
which momentum is transferred to the object, Q. The optical force, F, can be determined 
through 



F = Q 



(1) 
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where P is the trapping power, is the refractive index of the medium, and c is the speed 
of hght. We are interested in observing how force varies with position along the axis and as 
such output force curves that are either a function of lateral, or axial displacement. From 
these force curves several parameters can be taken or calculated that describe the system 



under study as re-illustrated in fig. IT] |28|. 



An object placed at a point where it experiences zero total force and is surrounded (within 
a certain proximity) by a negative gradient is said to be in equilibrium. Should the object 
be displaced the local gradient will produce a restoring force back toward this equilibrium 



position. In figures 1(a) and 1(b) these positions, z^q and Peq, are the axial and lateral 
equilibrium positions respectively. 

To achieve a single beam gradient force trap the efficiency, Q, in the negative axial 
direction must overcome that opposing it. If accomplished the efficiency curve will at some 
position become negative, allowing a point of zero force and negative gradient to exist. The 



maximum magnitude of this, Q^max^ is a good measure of the axial optical trap strength (28| . 
with its sign indicating whether a stable equilibrium position exists or not. 

The axial and lateral trap stiffness is proportional to the gradient of the force curves at 
their equilibrium position, which for the lateral case, assuming symmetry, is at zero. The 
trap stiffness, k, is related to the efficiency, Q, by 

UmPdQ 

^ = -JT- 2) 

C OS 

where s is either z or p for axial and lateral respectively. 

During the results and discussion we shall look at the force curves alone and also results 
that are determined from many such curves where the parameters have been varied. 

II. THE PROBLEM 

Ideally we want to model forces on a microdroplet trapped near the focus of a beam given 
only the properties readily known. The geometry of our system is significantly different to 



those normally studied [29|, [sQl as an additional discrete step in refractive index is located 



between the microscope objective and trapped particle due to the method of trap loading 



27l |. The detailed system is illustrated in fig. |2] with a laser beam of wavelength A, waist 



w, incident upon the back aperture of an objective lens of focal length /, aperture p and 
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FIG. 1: (Color online) (a) Axial efficiency, Qz, as a function of axial displacement. The axial 
equilibrium position, Zeq, occurs where the curve crosses y = with negative gradient, Qj^^nax 
quantifies the traps axial strength, (b) Lateral efficiency, Qp, as a function of lateral position. The 
lateral equilibrium position, p^q, occurs where the curve crosses y = with negative gradient. For 
(a) and (b) the gradient of the curve at Zqq and peq are proportional to the axial and lateral trap 
stiffness, and Kp respectively. 



focussed to the diffraction limit at a converging angle 9o, normally quantified in terms of 
numerical aperture, NA = n^sin^o- 

Our system assumes matching of objective lens, oil and glass coverslip refractive indices 
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(b) 

FIG. 2: Optical system and its parameters, a) A beam of wasit w enters an objective lens of focal 
length / with a back aperture of radius p. It is focussed to a point fp having propagated through 
two mismatched refractive index interfaces, zi and Z2, such that the thickness of the middle medium 
is Ah. If neither interface existed then the light would be focussed to point rf. b) Expanded view 
of the focal region of the microscope objective to the left. Light is incident on the first interface, 
zi, at an angle 9g and refracted to an angle 9w It is then incident on the second interface Z2, a 
distance Ah away, at an angle 6^ where it is refracted to an angle 6a and focussed to its parax;ial 
focus point fp. The height of the paraxial focus above the second and first interfaces is L2 and 
Li respectively. The droplet is trapped a distance h above the first interface, z above the paraxial 
focus and Zf below the point Vf. The first interface is between glass and water, and the second 
interface is between water and air. 



resulting in three consecutive layers of media and two refractive index mis-matched inter- 
faces. 



For the problem stated the decision of which size regime in which to calculate optical 
forces is important. For the particle radii here, i? ^ A, we are outside the Rayleigh regime 
and to a first approximation a GO approach would seem appropriate. However, considering 
the problem we see that this is not the case for the following reasons. 

Firstly, the NA, hence opening angle of optical traps is large and the paraxial as- 
sumption (sin^ ^ 9) is no longer applicable. For highly convergent beams the focus is 
not Gaussian but rather governed by an ir iteg ral representation due to the electromagnetic 



diffraction within the optical system [3l|, |32|. Also, as fig. |2] shows the beam is focussed 



through a coverslip and an aqueous layer. The interfaces created by this glass slide to aque- 
ous layer and aqueous layer to trapping medium (air) creates a mismatch in refractive index 
through which the beam is focussed. These discontinuities introduce spherical aberration 
into the focussed beam and so can only be fully described using a full wave analysis. 

Secondly, the interaction of a plane wave with a sphere where i? ~ A is more challenging 
to describe than by constructing the problem as a beam of many single rays passing through 
a sphere, as in GO. The description really must take into account diffraction. There is an 
analytical solution to Maxwell's equations for the scattering of a plane light wave by a single 
sphere for any ratio of radius to wavelength. The solution was independently developed 
by Mie, Debye, and Lorenz around the turn of the 20th century but has been historically 
referred to as 'Mie theory' or 'Mie scattering' {sSj. 

Thirdly, there is something intuitively wrong with the wavelength independence of GO. 
The focal waist and scattering of light by colloidal particles is known experimentally to be 
wavelength dependent so a theory of optical tweezers should also be wavelength dependent. 
An additional complexity arises from the proximity of reflecting surfaces which can cause 
'reverberations' [19] of light that significantly alter the interaction. 

Finally, a true description of the physics at play must traverse the full range of applica- 
bility from Rayleigh scattering to GO. 

III. THEORY 

We will initially utilise results from a GO approach to describe the forces on spheres to 



enable a comparison to a large number of relevant articles already in literature [34143 6l| . We 



chose to use the GO model, described in detail by Mazolli et al, which takes into account 
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the Gaussian nature of the beam and the Abbe sine condition 



37|. 



It will become apparent that theories beyond GO are needed to accurately model our 
problem. For this purpose we rely heavily on the Mie-Debye model given by Mazolli et 
al. S^] and its extension to Mie-Debye-Sphercal-Aberration (MDSA) theory described by 
Viana et al. [l^. Viana et a/.'s work includes a description of a beam after propagation 
through an interface of mismatching refractive indices. An extension of this theory to include 
two interfaces of mismatched refractive indices, accurately describing the geometry of our 
experiment, is now presented. 

We will first describe the tight focussing of a beam using the integral representation 
developed by Richards and Wolf 



31 



32| which can be extended to include propagation 



through stratified media 



38l |. Secondly we will give a brief overview of Mie's classic solution 
that describes light scattering from spheres and finally we will describe how optical forces 
can be calculated through combining these ideas. 



A. Focussed beam description 

A Gaussian laser beam with plane wavefronts entering the back aperture of a lens is 
described by 

^;obj(p,^) = i5obje*'"^e-^'/"'e, (3) 

where /cq = 27r/Ao, z is the axial direction, p is the lateral direction, and e is a unit vector 
along the wave propagation direction. 

When focussed into a medium of refractive index ng, in our case glass, the beam will 
occupy a conical region in space governed by the angle of convergence 60, and the azimuthal 
angle (p. This beam can be thought of as a superposition of plane waves and given by an 



integral representation of electromagnetic diffraction described by Richards and Wolf 
The electric field in glass is therefore 

f27r />eo 



3l| 



^^giass = ^o/ / sin^^,v/^e-^'^^^'^«e-^^-^^e*^e'(^,¥.)rf^t/¥;, (4) 
Jo Jo 



where 



71 f 

Eo = -^-f^Tobj^obj, (5) 

and 9o is the opening angle of the focussed beam given by the NA of a lens, 6*0 = 
arcsin (NA/ug), Ug is the refractive index of glass, 9g is the half-cone angle in glass, Tobj 



is the transmission of the objective, Eotj is the electric field magnitude at the objective lens 
back aperture, 7 = f/w and /, u and rj are defined in fig. [2l e{9,Lp) rotates the plane 



waves to occupy all angles from to 6*0 and all (p using a rotation by Euler angles [381, |39| 
{ip,6g, — v^). From here the NA of a lens describes the opening angle of the cone of focus 
exterior to and immediately before any interfaces. When the opening angle of the converging 
beam is larger than the critical angle for the glass to air refractive indices, the NA of the 
beam is effectively reduced to 60 = arcsin {N1N2). 

Equation 4 is the classic representation of a beam focussed to a point, Vf, however, our 
system differs as it has two interfaces between exit from the lens and reaching the focal 
point as shown in fig. [2J The plane wave components of equation H] each refract at the 
interfaces at zi = —h and Z2 = —h + Ah. From Snell's law the angle of refraction in the 
water layer 9^ = arcsin {sai9g/Ni) and the angle of refraction in air 6a = arcsin {sm6w/N2) 
where A^i = rising and N2 = ria/n^ are the relative refractive indices of the glass-water and 
water-air interfaces respectively. 

Our beam representation must therefore include the effects of propagation through media 



of stratified refractive index 

by 



. The focussed beam in the third medium, air, is described 



E^, = Eo r rTieg)smeg./^ge-^"''''"'^e-'^''^^-'-^''e-'''-^fe^ (6) 
Jo Jo 

where {kg^ — k^iz) h takes into account beam propagation in the glass slide up to the first 
interface, ka = riako, and each plane wave amplitude is multiphed by its respective Fresnel 
transmission coefficient 

TiOg) = UOgmOg) = ^^^^ eos 9^ COS el Tn2Cos 6, ' 

The effect of the additional factors in equation [6] over equation|4]is to introduce a spherical 

aberration that deforms the wavefront preventing diffraction limited focussing to the point 

Tf. This is quantified in terms of an aberration function as will be shown later in section [Tll Ci 
Having focussed the beam through two mismatched refractive index interfaces, the height 

of the paraxial focal plane above the water layer is found from the objective displacement, 

X, through 



(8) 



Using the work of Torok and Varga [38| and equation |6] we calculate the profiles of beams 
focussed in our system and compare them to the ideal beam assumed in most cases, giving 
some insight into the physics. 

Figure [3] displays the yz-plane beam profiles for focussing in water, through a glass- water 
interface, and through glass-aqueous-air interfaces. The axial displacement zero point is the 
position that the paraxial focus, Vf, would exist at when no refractive index interfaces are 
present. 




y (pm) y (nm) y (nm) 



(a) (b) (c) (d) 

FIG. 3: (Color online) Profile of a focussed 532 nm Gaussian beamstaken from a y-z slice through 
the beam axis, (a) The beam is focussed into water (n^ = 1.33). (b) The beam is focussed into 
water having crossed a glass (ug = 1.517) to water {riw = 1.33) interface after the lens, (c) The 
beam is focussed into air (ua = 1.00) across glass {ug = 1.517) to water (n^o = 1.342) and water 
to air interfaces. The objective displacement X = 40 ^m, the water layer is 10 /um thick, 7 = 1 
and ^0 = 41.23°. Zero on the axial axis is the position of the paraxial focus had there been no 
interfaces. 



The beam focussed in water with no preceding interfaces varies smoothly at the focus 
compared to those focussed in water and air having first travelled through glass coverslips. In 
particular the beam focussed to a point in air has a large number of oscillations in intensity 
along the beam axis. Previous work has shown such landscapes can interact with particles 
in a non-trivial manner |4oi]. Particles that are relatively large may not 'see' the oscillations 
while small particles could be trapped at more than one of the 'hot-spots'. The colour scale 
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remains the same for all plots so the maximum intensity is less in an airborne trap than for 
others given the same input power. 

Having described beam focussing more realistically, specifically for stratified media, we 
will move onto the scattering of light by the particles we wish to model in these beams. 



B. Mie scattering 



A full and clear derivation of Mie theory can be found in Bohren and Huffman 4l|] so 
here we only include the results that are important for our work. 

A plane wave incident on a spherical particle results in a scattered electric field that is 



dependent on the Mie scattering coefficients a„ and 6„ j4l| where n appears because of the 
Legendre polynomials in the solution which have n degrees. To simplify the following the 
Riccati-Bessel functions are introduced as; 

ipn{krsO = kvsJnikrs) and in{kr = kr sh^^\kr s) , (9) 

where k is the wavenumber, is from the spherical coordinates system and hn"^ = J„ + iyn 
is the spherical Hankel function with the spherical Bessel functions Jn and ?/„. Using these 
the scattering coefficients are 

and 



(10) 



where ^ = rimkoR is the size parameter and m = np/rim, the relative refractive index of 
particle to medium. 

The previous sections give the necessary background to now allow a description of how 
the force in optical traps is calculated. 



C. Force calculation 



To calculate the force, F, we follow the full electromagnetic approach, integrating the 
Maxwell stress tensor over the surface of the object; 
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h-TdS 



(12) 



where n is the surface normal and T is the electromagnetic stress tensor. Due to system 
symmetry and also momentum conservation the force simplifies so the surface is at infinity, 
thus giving 



(13) 



where E = Emc + Egcat, with equivalents for the magnetic field |37] . 

An analytical solution to this can be complicated, even for spheres, as seen in Barton 
et al. 42|. We take the simpler approach described by MazoUi et al. |37j where the vector 
electric and magnetic fields are given in terms of scalar Debye potentials, also known as Hertz 
vectors 43|. The optical forces are calculated by following Farsund and Felderhof 4^, who 
derive force, torque and absorbed energy for an object of arbitrary shape and material 
given the Debye potentials for the incident and scattered fields. The analytical solutions for 
force found through Farsund and Felderhof 4J] are converted to trapping efficiency through 



equation [Tl 

Debye 45|] calculated the force on a sphere due to an incident plane wave, and here the 
result is generalised to a focussed beam. First, the Debye potential for a single plane wave 
is 



/47r(2j + 1) 



j u 



E e-^^-^~'^^''dl^,{e,)Y,UO,v). (14) 



where Jj are the spherical Bessel functions and Yj^ are the spherical harmonics. Using the 
matrix elements of finite rotations, also known as Wigner d functions , for rotation 

in the basis of spherical harmonics, the Debye potential for a focussed Gaussian beam made 
from a superposition of plane waves, whose field is represented by equation [6l is 



6o 



nL(r,^,</^) = ^ / sin^fcv/^e-^'^-''^'= Ez^-V,(A;r 



X 



^ <i(4)y,„(0,^) Te-^-fe- 

Jo 



J U + 1) 



{kgz—kwz)h ^—i(m-l)ipk 



difk- (15) 



M=-j 
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Evaluating the integral over the azimuthal angle the Debye potential for the incident field 
becomes 



X 



m=-j 



X Jm-i{kpRsm9k)e 



From fig. [2] the relative locations of the planes gives 

Zf = — (Ah+^ 
^ nA N2 



Li — z. 



(16) 



(17) 



Substituting this into the middle two exponents of equation [16] we derive the aberration 



function 



19 



38 



39 



48j, \E', of our system to be 



^ = fco I - 



j^Ah + jfjfL2 1 cos 6g + n.^/\h cos Oy, + (L2 + z) cos 6'„ j . 



(18) 



The Debye potential for the scattered field is found through the incident fields interaction 
with a sphere and hence is dependent on Mie coefficient and Hankel function h^p such 
that 

00 j 



im-l)(pf 



j=l m=-j 



IAtt (2j + 1) 

J U + 1) 

xa,hf\kr)Y,m{0,v), (19) 



where 



Gjm= T{9) sineV^e-^ '''' '^di^M-Jm-iikpsmea)e'^^'^'^^de. (20) 
Jo 

Similar expressions can be found for the magnetic field, H, using the Mie coefficient 

The efficiencies are given for the lateral and axial components each with two separate 
contributions, one for the rate of removal of momentum from the incident beam, Qe, and 
the other for minus the rate of momentum transfer to the scattered field, Qg, so the total 
efficiency Q^^^ = Q^'^ + Q^'^. The forces are calculated for circularly polarised light but can 
equally but done for linear polarisations 48|]. The axial component of the trapping efficiency 



is given by 



19 



37 



48 



q: 



47^ 

AN1N2 



00 J 



'^Y.Y. (2j + l)(a, + &,)G,-^G;* 



(21) 



j=l m=-j 
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and 



q: 



j=l m=-j 



J + 1 



(22) 



J (j + 1) 



The lateral efficiencies are 

2 oo J 



27^ 



" AN1N2 



(23) 



and 



Q 



' V- (j + 2) (j - m + 1) (j + m + 1) , ^ ^ ,^ . 

j=l m=-j 

X (^i,mG'j+l,.m+l + ^j,-m^i+l,-m-l) ' (24) 



where A is the fraction of the beam power that enters the objective back aperture, to account 
for overfilling 19[, and the functions G'j „^ and are defined as 



G]^ = j\{e) sin0v^cos^,e-^'^^'^'V^,i(0JJ„_i(fcpsin0Je^*(^'^)rf^. 



(25) 



and 



G 



jm 



/ T{e) sin sin ^.e"^' '^dl^^^^ (^a)^m-i(A:p sin ^Je^^^^'^'^d^. (26) 

Jo 



In the limiting case where Ah = and riy^ = Ug the results return to those of Viana et 
al. 1^ for a glass to water interface without an intermediate aqueous layer. For A/i = and 



ng = Ua = Uw and X = 0, the results of Mazolli et al. 



371], and Neto and Nussenzveig [4 



are matched. 

The final crucial point concerns the computation of these equations. Rather than com- 
pleting the sums in equations [21] to |24] ad infinitum it is useful to know that it is sufficient 

1/3 

to sum over j up to S,s + + 2, or its nearest integer, due to the convergence of the Mie 
scattering coefficients (Appendix A in Bohren and Huffman 41|). 

Having introduced the theory and visualised the focussed beams we will move onto ex- 
amining the outcome of applying the theories to 'normal' optical traps and then to our 
application. 
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IV. RESULTS AND DISCUSSION 



Ashkin observed that even with relatively loose focussing of a Gaussian beam, particles 
(with m > 1) always had the tendency to move toward the beam centre where they would 
reach a lateral equilibrium position. In airborne tweezing this occurs in exactly the same 
manner so its modelling is not of great importance. It is the axial efficiency and force curves, 
and the associated balance between gradient and scattering forces, that governs whether a 
particle is trapped or not. Although for completeness we have given the analytical solutions 
for lateral efficiencies, it will be the axial direction we consider as it determines the unusual 
phenomena observed. 

The results we present should be considered a set of typical examples that can be produced 
using our code and is by no means exhaustive. We also note that we have included several 
results pertaining to the optical trapping of spheres in water to highlight the large difference 
between those experiments and those in air. 

A. Comparison of geometrical optics and Mie scattering 

First we will make a comparison between the theoretical predictions of GO against those 
from Mie scattering. In this first instance we will neglect the effects of spherical aberration 
and show in fig. H] the axial trapping efficiency calculated through both theories when a 
250 nm, 1 fim, and 5 fim silica sphere is trapped with 532 nm light in water. As described 
in the theory section there is a limit to the opening angle of the focussed light and hence NA 
of the trapping beam. For a beam focussed through a coverslip-water interface this limit is 
6o = 6c — 62°, thus NA = 1.33 sin = 1-17, and is the value used for this first test. 

Clearly the prediction of GO disagrees with those of Mie scattering. However, GO stands 
up surprisingly well even for spheres with radii similar to the wavelength of trapping light. 
Testing the theory on a 5 fim sphere, which is approaching the regime where GO should 
become applicable, it is indeed a reasonable approximation except for the area closest to the 
paraxial focus and at the extremities. 

The inaccuracies are unsurprising considering the earlier discussion of the limits of GO. 
For small spheres Mie scattering plays a dominant role that differs to simple ray optics and 
for the larger sphere the non-Gaussian beam focus plays the important role that GO cannot 
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FIG. 4: (Color online) Axial trapping efficiency calculated through GO and Mie theories when 
250 nm, 1 /im, and 5 /im silica spheres {up = 1.445) are trapped with 532 nm light in a water 
medium (n^ = 1.33) with 7 = 1 and = 61.25° in a system like fig. [2] with no refractive index 
interfaces. The four curves are plotted on two separate graphs for clarity. In (a) the black solid 
line is calculated through GO and the blue dashed line is calculated through Mie scattering. In 
(b) all curves are calculated with Mie scattering. 

account for. 

We will now test how applicable GO is when trapping objects in air. Here the upper 
limit of the NA is reduced (^0 = — 41.2° therefore NA ~ 0.66) and lends itself toward the 
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paraxial approximation, hence GO. Yet the ratio of particle to medium refractive index is 
higher than in colloidal systems thus moving further into the applicability of Mie scattering. 
Figure plots the axial efficiency for the same particles as fig. H] except the medium is now 
air [rim = 1.00) and the particle is a water droplet {up = 1.342). 

The form of the GO curve remains similar to the previous example except with an absolute 
increase in efficiency. This is probably due to an increased scattering force component from 
the larger particle-medium refractive index contrast and reduction in NA. GO predicts the 
droplet will 'just' not obtain an axial equilibrium position allowing it a brief reprieve in 
matching the more rigorous Mie theory. However, this disappears quite quickly when noting 
the drastic curve change as the forces on three sizes of spheres are computed using Mie 
scattering. The largest sphere, 5 fim, enters the beginning of the GO regime {R^ A), yet 
the theory completely fails to indicate the occurrence of a second minima, predicted by Mie 
theory. 

It has been shown here that GO, although not highly accurate, can provide reasonable 
predictions of the efficiencies of trapping colloidal particles in water, giving indications of 
what one would expect in real world systems. However, in the same manner it has been 
shown that GO is not an appropriate description of airborne trapping with the wild variation 
as a function of radius not predicted, and the inability to predict important features. 

Our next extension is where the GO description falls down, namely in the consideration 
of the relevance of spherical aberration for both colloidal and airborne systems. Only Mie 
theory can accommodate variations in the phase wavefront profile so now only this will be 
considered. 

Figure [6] plots the axial efficiency curves for 1 /im and 5 /im silica spheres trapped in a 
water medium when the aberration due to a single coverslip-water interface is and is not 
neglected. 

Spherical aberration clearly has a significant effect on optical trapping efficiency curves as 
expected 50, ImI]- There is a drastic difference between the axial efficiencies from aberrated 
and non-aberrated beams. The two main effects are a reduction in Q^max^ reducing the axial 
trap strength, and a decrease in axial equilibrium position, Zeq, so the spheres 'sit' lower in 
the trap relative to the paraxial focal point. 

Spherical aberration plays a major role in the physics describing optical trapping, so it 
must be considered in our system (fig. [2]) where there are two interfaces with mismatched 
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FIG. 5: (Color online) Axial trapping efficiency calculated through GO and Mie theories for 250 nm, 
1 fim, and 5 fim water droplets (n^ = 1.342) trapped with 532 nm light in air (rim = 1.000) with 
7 = 1 and 6q = 41.23° in a system like fig. [2] with no refractive index interfaces. The four curves 
are plotted on two separate graphs for clarity. In (a) the black solid line is calculated through GO 
and the purple dashed line is calculated through Mie theory. In (b) both curves are calculated 
with Mie theory. 

refractive indices. Repeating the previous figure for water droplets trapped in air above a 
coverslip and thin aqueous layer we obtain fig. [71 

The inclusion of spherical aberration in the description greatly affects the efficiency curves 
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FIG. 6: (Color online) Axial efficiency curves for 1 ^um and 5 /xm silica spheres {up = 1.445) 
trapped in water above a glass coverslip with and without the aberration induced by the refractive 
index interface of glass {ug = 1.517) to water (n^ = 1.33) taken into account. The objective axial 
displacement X = 35 /um, 7 = 1 and 9q = 61.25°. (a) For a 1 /um sphere the blue solid line is 
without aberration and the purple dashed line with aberration, (b) For a 5 fim sphere the red solid 
line is without aberration and the black dashed line with aberration. 

for airborne water droplets. There is a reduction in Q'^axy reducing the axial trap strength 
and for the larger droplet a general 'smoothing' of the curve occurs with smaller local minima 
created. Figure El for the same system and objective displacement but with and without the 
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FIG. 7: (Color online) Axial efficiency curves for 1 /xm and 5 /im water droplets (n^ = 1.342) 
trapped in air (ua = 1.000) above a glass coverslip {ug = 1.517) and thin water layer (n^ = 1.342) 
as depicted in fig. [2j The objective axial displacement X = 40 /im, the water layer is 10 fim thick, 
7 = 1 and = 41.23°. (a) For a 1 /um sphere the blue solid line is without aberration and the 
purple dashed line with aberration, (b) For a 5 iim sphere the red solid line is without aberration 
and the black dashed line with aberration. 

aqueous layer, shows a change in the axial curves indicating the significance of the aqueous 
layer and its necessary inclusion in the theory. 

Having established the most appropriate model to use and the physics to include - pri- 
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FIG. 8: (Color online) Axial efficiency curves for 1 /xm and 5 /um water droplets (n^ = 1.342) 
trapped in air (n^ = 1.000) above a glass coverslip {ug = 1.517) with and without a thin water 
layer (n^ = 1.342). The objective axial displacement, X = 40 /im, 7 = 1, = 41.23° and when 
the thin water layer exists it is 10 /xm thick. In (a) the blue solid and purple dashed curves are 
calculated without and with the thin water layer respectively. In (b) the red solid and black dashed 
curves are calculated without and with the thin water layer respectively. 

marily the effect of the stratified layers between the objective and trapped particle, spherical 
aberration and the high relative refractive index - we will now move on to attempt to explain 
the phenomena observed as described in the introduction. 
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B. Predicting experimental observations 



All the previous theoretical results shown have neglected any particle buoyancy. For 
colloidal systems this is a reasonable approximation with the density of the trapped objects 
approximately that of the medium. Thus, there is only a multiplicative factor between 
efficiency and force graphs via equation [1] which allows the axial Q curves to be treated 
as scaled force curves. However, this is a very poor assumption when considering water 
droplets suspended in air with the large density contrast: Pwater — lOOOpair- In order to 
fully appreciate what the theory predicts we must calculate the force experienced by the 
microsphere using equation [T] and subtract its weight. The droplet and system parameters 
from fig. [7] including aberration is replicated with a trapping power of 10 mW to calculate 
the corresponding force curves in fig. [91 




FIG. 9: (Color online) Axial force curves for 1 ^um and 5 /im water droplets {up = 1.342) trapped 
in air [ua = 1.000) above a glass coverslip {ug = 1.517) and thin water layer (n^ = 1.342). The 
objective axial displacement, X = 40 /im, the water layer is 10 ^m. thick, 7 = 1, = 41.23°. The 
red dashed line is for a 1 /xm droplet and the blue solid line is for a 5 ^m. droplet. 



For the 1 pm, unlike the 5 /im, droplet the force plot has no significant effect on the 
properties deduced from the efficiency plot. The weight of the 5 /im droplet is comparable 
to the optical force so that an equilibrium position, Zeq., exists which was not indicated in 
the optical efficiency curve. 

We now examine if the theory predicts some of the behaviour observed during aerosol 
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trapping experiments starting with points one and two in the introduction. In fig. 10 (a' 



we 



plot for a 4 fim water droplet, trapped in the experimental system depicted in fig. |2l the 
predicted axial force curves for increasing trapping powers. Repeating for several droplet 
radii, the height above the water layer that a droplet is trapped, obtained from Zf,q, can be 



plotted as a function of power as shown in fig. 10(b) 



Figure 10(a) successfully predicts two physical observations from experiments. As the 
trapping power increases the droplet equilibrium position, Zeq, and hence height above the 
underlying water layer, increases and with enough power eventually falls from the trap 22 |. 



The curves of 2 /xm, 3 /xm and 4 /xm droplets in fig. 10(b) do not continue for higher powers 
but the 5.5 fim droplet continues indefinitely. Clearly, if an equilibrium position exists in 
the efficiency curves alone, then the droplet will always remain trapped. If no such position 
exists then the force curve may eventually lose its equilibrium position with increasing power. 
This qualitatively explains our own experimentally observed results and the power gradients 



of Knox et al. 



23|, indicating their measured gradients [23 1 gment of the extended 



curves. With this enhanced understanding their suggestion of using power gradients as a 
tool for aerosol sizing could benefit. 



Figure 10(b) may also explain point four, why there is a linear dependence on captured 
droplet size as a function of power, and why small droplets cannot be trapped at high powers 
(figure 4 in ^| and figure 2 in [24] )• The power gradients show that above certain powers, 
depending on droplet radius, it is possible for no equilibrium position to exist. Therefore, 
although a 'large' droplet may be trapped at relatively large powers, smaller droplets cannot 
be for the same power. 

A large parameter that governs the magnitude of the spherical aberration induced by the 
interfaces is the depth into the sample which the beam is focussed. For example, a lower 
focus has less aberration. In fig. [11] the beam is simulated to focus at several depths into 
the sample chamber and the force curve calculated again for a 4 /im water droplet. 

The decrease in aberration not only shifts Zeq closer to the paraxial focus but also increases 
the strength of the optical trap with an increasing F~^^^ and overall deepening of the 
potential well. 

Now consider point three from the introduction. It is reasonable to expect the particle to 
exhibit some sort of interferometric properties with refiections from the inside of the droplet 
interfering with themselves. This can be demonstrated by simplifying the model and treating 
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it as an interferometer which performs reasonably well at estimating the axial efficiency at 
the paraxial focus 49|]. Liquid aerosols will establish a stable size once in equilibrium with 
their surrounding environment, namely the relative humidity [81]. Although the process of 
growth and evaporation is relatively fast it is at times clear one of these occurs just after 
the droplet becomes trapped. Investigating how the height at which the droplet is trapped 
varies with droplet radius we plot fig. [T2l 

There is a clear, near sinusoidal, oscillation in droplet height as a function of its radius. 
A single oscillation in height occurs over a change in droplet radius of ~ 100 nm, going from 
a local minima to maxima in half this, ~ 50 nm. So, a change in trap height of ~ 2 /xm 
occurs due to only a 50 nm change in droplet radius. When observing a particle just after 
capture the change in size is clear, far above the limit of resolution, so must be greater than 
50 nm. Knowing that the oscillations are most frequent just after capture we conclude the 
multiple oscillations that occur in experiments are due to changing particle radius and hence 
equihbrium position, Zeq- 

To measure this oscillation would be a challenging experiment. The droplets would need 
to be imaged from the side to measure their height and also coupled with a high precision 
sizing technique such as CERS [ref]. With such instruments in place the droplet radius would 
have to be varied by altering ambient relative humidity or varying droplet temperature [ref]. 
Both would also change the salt concentration of the droplet hence refractive index so this 
effect would need to be included in the model. 



C. Limits of techniques 

Section B qualitatively explained the appearance of four unique phenomena observed 
when trapping airborne water droplets by simulating an experimental system through mod- 
elling the optical forces created by the focussing of a high NA beam through two refractive 
index mismatched interfaces. In what follows we shall see how far the boundaries of optical 
trapping in air can be pushed. We shall investigate if smaller particles (< 1 /im) can be 
trapped, if the axial trap strength and capture volume can be increased, and determine 
the limits on the particle refractive index that can be trapped. To explore these questions 



Qz max is calculated as a function of both particle radius and relative refractive index [28| , 



first for particles suspended in water, fig. [131 ^ means of comparison, and then for airborne 
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particles; fig. [T] 

The white areas on the plots represent parameter space where a negative Qaxiai value 
does not exist and hence no stable trap position is possible 56|]. Of immediate note are the 
'spikes' in the contour plots indicating resonances in the force experienced by the particles. 
The effect is more pronounced as a function of radius although at the high refractive index 
end of the spikes there are rapid resonances in force as a function of refractive index, creating 
tiny islands of parameter space where traps can exist. 

These resonances can be explained by interference effects due to the spheres increased 



refelctivity at high relative refractive index and its associated variation with radius [28| . The 
decreased frequency of the resonances in air, fig. (TU is due to the lower medium refractive 
index (ra^ = 1.00 and = 1.33). 

As the plots are functions of relative refractive index it is noted that for a given particle 
refractive index the horizontal line of interest is higher up the refractive index axis in air 
than water. Looking at the sorts of particles normally trapped in both media gives a 'feel' 
for the plots. For example, a silica sphere in water exists along the line defined by rirei — 1-09 
in fig. [13] and for a water droplet in air the line is at Urei — 1-34 in fig. [TU 

The continuous region of stability for optical tweezers in air is over a smaller range of 
refractive indices {Arirei — 1 — 1.25) than when trapping in water (An^e/ — 1.33 — 1.65) and, 
also, the maximum negative axial efficiency values, Q^rnax^ smaller overall for trapping in 
air than in water. This is understandable because the larger relative refractive index between 
particle and medium for aerosols increases the Fresnel reflection coefficients, hence increases 
scattering forces which probably overcome the gradient forces. The minimum radius possible 
to trap is smaller in water than air probably due to the increased spherical aberration in 
the focussed beam, induced by the coverslip interface, which has a larger refractive index 
contrast in airborne traps. As seen in fig. [3] this increased aberration produces larger period 
oscillations in intensity allowing more 'room' for particles to 'fall' between. 

These plots are deceptive; the true range of particles that would theoretically obtain an 
axial equilibrium position in air has been misquoted because, as previously mentioned, the 
weight of the particle is significant. F~^^^ is the truly relevant quantity that will allow the 
determination of whether the spheres are isolated in three dimensions or not. Unfortunately, 
this poses a problem as the force from equation [1] is dependent on laser power and with this 
additional variable not all parameter space can be easily displayed. Instead F~^^^ is plotted 
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for a single power, P = 10 mW, in fig. [151 

Comparing figures [H] and [15] we come to an interesting conclusion. For droplets with 
certain particle parameters, indicated in fig. [TU traps are created through the transfer of 
optical momentum alone (single beam gradient force trap or optical tweezers). However, 
fig. [15] indicates that with the assistance of gravity a larger range of droplets can be 'trapped', 
although not with momentum transfer alone. Consider a droplet that evolves in size (it will 
also evolve slightly in refractive index due to salt concentration changes); as the radius varies 
the 'path' of the particle in the parameter space of fig. [H] may cross through a non-tweezing 
region but due to its weight remains trapped (fig. [T^ . This difficulty in deciding whether 
a droplet is tweezed or levitated leads to the conclusion that as a general name for the 
experimental tool being used we really have a quasi optical tweezers. 

In fig. [16] we superimpose the tweezing and trapping areas of figures [H] and [15] Areas of 
parameter space truly optically tweezed are coloured grey, areas that are only trapped with 
the assistance of gravity are coloured red and the area that would be truly tweezed if the 
droplets had neutral buoyancy is coloured blue. White areas retain the same meaning of 
neither tweezing nor levitation. 

It is clear now that the choice of inverted or non-inverted tweezers is critical in the success 
of optically trapping a large range of aerosol sizes. Having established the true nature of the 
technique we are using the next section discusses the possibility of moving more of parameter 
space into the tweezing regime (trapping through momentum transfer alone). 

D. Optimisation and extension of limits 

We have demonstrated several points that stop airborne traps from reaching their opti- 
mum performance. These include spherical aberration created in the beam, a high refractive 
index contrast between particle and medium causing large scattering forces, and a lack of 
high converging angles (reduced NA) due to total internal refiection at the coverslip interface. 

Total internal refiection is not easily circumvented but it could be possible to correct 
for spherical aberration or possibly remove the large scattering forces which we shall now 
explore. 
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1. Spherical aberration correction 



It is feasible to correct for the spherical aberration using spatial light modulators, de- 
formable mirrors, or other wavefront modfiying elements {52 1. Any correction would clearly 
be advantageous creating better localisation of aerosols and hopefully moving into the im- 
portant accumulation mode size regime. In fig. [17] we plot for the same parameters as fig. [H] 
but an additional spherical aberration, of magnitude 0.08A, is placed on the beam input to 
the objective. 

Figure [T71 shows an improvement in the axial strength of the optical tweezers, an increase 
in the overall range of parameter space that can be tweezed but unfortunately also an 
increase in the minimum particle radius that can be tweezed. 



2. Removal of beam centre intensity 

Ashkin used GO to predict an increase in Q~^a.x if ^ 'doughnut' mode beam (TEMqi) 
fills the back aperture of an objective. Increasing Q^max is by far the most difficult problem 
in airborne tweezers, as shown and discussed, so next we predict the effects of removing 
the central portion of a Gaussian beam as it enters the objective. In fig. [18] we plot Q^max 
against both radius and relative refractive index for a Gaussian beam where ~ 57% of the 
beam area is removed leaving an annulus, yet keeping the total power the same. 

Figure [TS] shows that the area in parameter space over which a true optical tweezer can 
be created is greatly increased by removing the central core of a Gaussian beam, although 
the minimum sphere radius tweezable has increased. The minimum radius increase could be 
due to the zero intensity that may now exist in the focal plane of the tweezers into which a 
small enough sphere could sit, experiencing no forces from the surrounding light. 

The predicted increase in parameter space over which aerosols can be tweezed is of great 
promise to the field of aerosol optical manipulation. It is difficult to trap high refractive 
index aerosols, specifically solid microspheres, yet they are of great importance to the fields 
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15| . It is hoped that a 



of medicinal drug lung delivery and atmospheric chemistry 
definitive experiment can be performed in future to verify this increase in optical tweezers 
parameter space. 
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E. Capture volume 



We have demonstrated that power gradients may explain the hnear dependence of cap- 
tured droplet size on trapping power but there may be more interactions occurring than 
thought. In order for the droplets to become trapped they must enter a capture volume, so 
it would be pertinent to calculate how this volume varies with trapping power and droplet 
radius. This volume extends between the maximum and minimum force points in the ax- 
ial and lateral directions simultaneously. Unfortunately, evaluation of this volume requires 
calculation of forces for locations away from the optical axis. At these locations there is 
a complex interplay between axial and lateral efficiencies [lol, [s^]. This will require more 
study to ascertain a suitable description and answer to the question of whether the capture 
volume plays a significant role in the linear dependence of captured droplet size on trapping 
power. 

V. SHORTCOMINGS OF THEORY 

Within the Mie scattering theory outlined above we make the assumption that the Fresnel 
transmission coefficients, tg and tp, for TE and TM modes of polarisation respectively are 
equal and take tg = tp. Plotted in fig. [19] are both modes of the transmission coefficients 
as a function of incident angle on a glass to air interface to highlight the inaccuracy in this 
assumption. 

The discrepancy in fig. [19] is significant. We propose that in future work it should be 
attempted to place both TE and TM modes into the theory. 

VI. CONCLUSION 

We have outlined a modified MDSA model of optical traps to describe our experimental 
system, and shown that it must be used over a simpler GO approximation in order to fully 
describe the beam profile and scattering. Investigation of the MDSA parameters indicate 
spherical aberration must be included, and we have discussed that gravity must not be 
naively ignored. 

Once correctly described we use our extension to MDSA theory to qualitatively predict 
and explain four unusual phenomena observed only in airborne optical traps. It has also 
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allowed us to explore limitations of current experiments and approaches before investigating 
what can be improved. 

The resonance type plots go some way to explaining what the experimenter observes: 
it has been observed experimentally that water droplets [rirei — 1.342) are easily trapped 
for a wide range of sizes, but we can say from experience that it is difficult to trap small 
water aerosols (< 1 /im) even though they are produced from the nebuliser. In our previous 
work js^l where solid aerosols [rirei — 1.445) are trapped, we can say that the expected range 
of particles trapped was odd. In colloidal systems, if two particle sizes can be trapped, almost 
certainly a size between these two will also be trapped. However, in air, spheres with certain 
radii could not be trapped yet sizes both above and below could. 

Looking at fig. [T5]both these phenomena can be qualitatively explained with the existence 
of the resonances as a function of radius and the lack of Q'^max small spheres (< 1 fim). 
Obviously, it is not easy to prove the non-result of being unable to trap certain objects, but 
the results here give some indication as to why it is so hard to trap small spheres, with a 
relatively high refractive index, that are so easily trapped in water. 

This work has lead to many new insights into how aerosols are trapped in single beam 
gradient force traps. Namely that an optical trap in the inverted geometry behaves as a 
quasi optical tweezer, at times tweezing droplets and at others only levitating them. We have 
given qualitative predictions that explain physical phenomena observed experimentally and 
explored the theoretical limits of trapping aerosols, both which help to define the parameters 
of the current tools at our disposal. The challenge for the future is to produce quantitative 
agreement between experiment and theory. 
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FIG. 10: (Color online) (a) Variation of axial force for a 4 /xm water droplet (rip = 1.342) trapped 
in air (tt-q = 1.000) at trapping powers of 10 mW (blue dot-dashed line), 20 mW (green dotted 
line), 50 mW (purple dashed line) and 200 mW (red solid line). The force for each power has been 
normalised to unity for clarity, as it is only the axial equilibrium position, Zgg, that is of conern. 
(b) Variation in height above the water layer droplets of radius 2 (blue dotted line), 3 /im (red 
dot-dashed line), 4 /um (purple dashed line) and 5.5 yum (black solid line) are trapped as a function 
of trapping power (power gradients of Knox et al. [23]). All but the 5.5 /um droplet curve stop due 
to the loss of axial equilibrium position at high powers as in (a). The objective axial displacement, 
X = 40 fim, the water layer (uy^ = 1.342) is 10 fim thick, j = 1, 9q = 41.23° and the coverslip 
refractive index Ug = 1.517 for both (a) and (b). 
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FIG. 11: (Color online) Variation of axial force for a 4 /xm water droplet {up = 1.342) trapped in 
air (ria = 1.000) with 8 mW of power for microscope objective displacements of 25 fim (red dot- 
dashed line), 30 ^m (solid black line) and 35 ^m (dashed blue line). The water layer (n^ = 1.342) 
is 10 //m thick, 7 = 1, 0o = 41.23° and the coverslip refractive index ng = 1.517. 




FIG. 12: (Color online) Plot of the height a water droplet {up = 1.342) in air (n^ = 1.000) is 
trapped above the underlying water layer (n^ = 1.342) as a function of radius. The objective 
axial displacement, X = 25 fj,m, the water layer is 10 /um thick, 7 = 1, = 41.23°, the coverslip 
refractive index Ug = 1.517 and the trap power is 10 mW 
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FIG. 13: (Color online)Q~^Qj. as a function of relative refractive index and radius for spheres 
trapped in a water medium (n^ = 1.33). The objective axial displacement, X = 40 /im, 7 = 1, 
^0 = 61.25° and the coverslip refractive index Ug = 1.517. The colour bar in (b) is representative 
for both plots, (a) is an expanded view of the first 2.5 /xm of (b) as this is the region where the 
boundaries would really like to be pushed. 
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FIG. 14: (Color online) Qj^^j^ as a function of relative refractive index and radius for spheres 
trapped in an air medium (n^ = 1.000). The objective axial displacement, X = 40 fim, the water 
layer (n^^, = 1.342) is 10 /im thick, j = 1, 9q = 41.23° and the coverslip refractive index ng = 1.517. 
The colour bar in (b) is representative for both plots, (a) Is an expanded view of the first 2.5 /um 
of (b) as this is the region where the boundaries would really like to be pushed. 
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FIG. 15: (Color online)^"^^^, as a function of relative refractive index and radius for spheres 
trapped in an air medium (ua = 1.000). The objective axial displacement, X = 40 fim, the water 
layer (uy^ = 1.342) is 10 fiui thick, 7 = 1, = 41.23° and the coverslip refractive index Ug = 1.517. 
The colour bar in (b) is representative for both plots, (a) Is an expanded view of the first 2.5 /um 
of (b). 
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FIG. 16: (Color online) Superposition of figures [T3] and [15] highlighting the areas of parameter 
space, as a function of relative refractive index and radius, where water droplets are truly optically 
tweezed (high R, low n^e/, grey), only trapped with the assistance of gravity (high R, high Urei, 
red), and optically tweezed if the droplet had neutral buoyancy (low R, low n.rei, blue). The white 
area represents areas where neither optical tweezing nor levitation occurs. The parameters for 
these plots are the same as the respective figures. 
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FIG. 17: (Color online) Q^^nax ^ ^ function of relative refractive index and radius for spheres 
trapped in an air medium (ua = 1.00) with a Gaussian beam entering the back aperture of the 
microscope objective with an additional spherical aberration placed on the beam at the entrance 
to the objective back aperture of magnitude 0.08A. The objective axial displacement, X = 40 /xm, 
the water layer is 10 fj,m (nyj = 1.342) thick, 7 = 1, = 41.23° and the coverslip refractive index 
Ug = 1.517. (a) is an expanded view of the first 2.5 fini of (b). 
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FIG. 18: (Color online) Qj^ax ^ function of relative refractive index and radius for spheres 
trapped in an air medium (ua = 1.000) with a Gaussian beam entering the back aperture of the 
microscope objective with 57% of its central area removed. The objective axial displacement, 
X = 40 fim, the water layer {riui = 1.342) is 10 /um thick, j = 1, 9q = 41.23° and the coverslip 
refractive index Ug = 1.517. 




FIG. 19: (Color online) Fresnel transmission coefficients tg (black solid) and tp (red dashed) for 
TE and TM modes of polarisation respectively as a function of incident angle 9 up to the critical 
angle ~ 41.5° for a glass (n^ = 1.517) to air (n^ = 1.000) refractive index interface. 
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